A Minimal Growth Equation For Conservative Non Equilibrium Growth From Vapor 
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A continuum growth equation in 1+1 dimensions is obtained by considering the contributions to 
the underlying current from kinetics of adatoms on the surface . These considerations reproduce 
the slope dependent term that represents a stable or an unstable growth depending on the sign of 
the associated coefficient. In addition an /i — » ~h symmetry breaking nonlinear term is obtained. 
Implications of this term on the stabler unstable growth transition are discussed. This term is 
tilt independent and is expected to be present independent of the source of the current in the 
conserved growth, computer simulation based on a zero diffusion bias model is used to verify 
some of the predictions of the minimal growth equation. It is shown that the growth equations 
for cellular automata type models involving instantaneous relaxations can be obtained from the 
proposed method by addressing the corresponding kinetics of adatoms. 
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A continuum description of conserved growth from 
vapor is obtained from Langevin type equation dth + 
W ■ J — F where /i(r, t) is height function , J is current 
due to adatom relaxation on the growing surface and F is 
the average flux with white noise. The current J includes 
an equilibrium and a nonequilibrium(NE) contribution. 
Most of the experimental growth is far from equilibrium 
and has been studied over last many years . Descrip- 
tion of such growth is essentially through the noncquilib- 
rium current Jne- The contributing terms to 3ne are 
obtained by comparing the behaviour of growth equa- 
tion containing sets of different terms, satisfying required 
symmetry and invariances of the problem, with the ex- 
perimental observations or, appealing to the un- 
derlying kinetics (^J^- Former is an empirical approach 
while later, a direct approach. It is desired that from 
kinetic considerations , Jne could be constructed, de- 
scribing the experimental observations. Such an approach 
was adopted by Villain earlier. In reference from 
computer based solutions of the equations describing step 
velocities in 1-1-1 dimensions, and the consequent inter- 
face profiles, it was concluded that the local form of the 
growth equation should comprise of a) a term of the form 
where a' and b' are constants, related to 
-h symmetry breaking 



(l-|-a'|m|)(l+fc'|m|) 

the diffusion kinetics, b) an h 
term and c) a term of the form V^ft.. 

In what follows we show that these terms can be de- 
rived by appropriately identifying the kinetic processes 
that contribute to the nonequilibrium current. To our 
knowledge the method proposed in this work has not 
been used to obtain J7v_e in any of the previous work 
related to growth from vapor. Most of our discussion 
will be pivoted around the symmetry breaking term that 
assumes the form V^(Vft.)'^ for small inclinations while 
at large slopes takes the form +-9a;m thus unifying the 
descriptions [^,^ at small and large slopes. While ob- 
taining the NE contributions to the current, we will con- 
sider growth with Schwoebel-Ehrlich (SE) barrier to 
account the stable as well as unstable growth. However 



the method of obtaining these contributions can be ap- 
plied to any other models not employing the SE barrier 
. We have applied it later to a stochastic model with 
cellular automata type rules that mimic the molecular 
beam epitaxial (MBE) growth at low temperature |ic| ] 
referred in the literature as DT model. 

Consider growth on a 1 dimensional flat substrate with 
lattice constant a. We will consider the situation depicted 
in Fig. |l| for obtaining various contributions to the cur- 
rent where steps are such that positive slope is obtained. 
Adatoms are randomly deposited on the substrate. Let 
Dg be the diffusion constant on the terrace, Ic be the aver- 
age distance that an adatom travels on a terrace before 
encountering another adatom. Detachment from steps 
or nuclei on the terrace is negligible. Adatoms hopping 
down the descending steps contribute to the downward 
current jd while those hopping on the terrace can get at- 
tached to an ascending step contribute to the in-plane 
current ji. Referring to Fig. |l| , the adatoms reaching 
site A and hopping down the step to the left constitute 
jd and those reaching site B and hopping to the right to- 
wards the ascending step constitute ji. The net current 
is jd + ji. These are obtained as follows, 

Jd(i)=(local density of site A(B))-(flux of adatoms 
approaching A(B))- (probability for hopping across 
A(B)). 

(1) 

Since the relaxation of adatoms is through the diffu- 
sion on the terrace or across the step edges, the density 
of sites A and B is same as density of steps which is ap- 
proximately given by j^j^, where m is the local slope. 

In the absence of nucleation, lateral flux(LF) approach- 



ing site B or A is 



±nF 



being the average local 



2|m|a-i ' 

terrace width and n is unit vector in the x direction. 
When terrace size is large enough, this flux is restricted 
due to the nucleation. The nucleation process will restrict 
the diffusion to an average length Ic on a large terrace. 
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As a result, even for very large terraces, the approaching 
flux LF is almost constant. The effect of nucleation is 
incorporated by introducing Ic in the expression for LF 

as 1 "-^ , — 7T, so that for small slopes the expression 

Let Pa and Pb represent 



fiF 

2(i-^ + |m|a-i) 

is reduced to a constant value 



the probabilities of hopping across sites A and B respec- 
tively. The Schwoebel length [|| Is oc {Pb — Pa) ■ Hence 
from the Eq.(l) the current is 



h\m\F{PB - Pa) 
2{1 + \m,\){lc^ + \m\a'^) 



(2) 



The LF approaching sites A and B is however modified 
due to the relative motion of the steps. Consider the sit- 
uation in Fig.0 where v is the velocity of the step bearing 
the terrace while v' that of the higher one on the positive 
slope. For v' > v the terrace width is reduced, depleting 
the LF approaching sites A and B. The reduction in flux 
is (X 5v where 5v = (v' — v). Adatoms hopping across 
upper step as well as those attaching in-plane, both con- 
tribute to the velocity of the step. Thus the velocity 
V oc js except that the coefficient (Pb—Pa) is replaced by 
(Pa + Pb). Hence Sv oc ,if-"+-^^\ 9^ , , , 

^ ' l^'- + \m\a-^ 2(l+|m|)(/^ H|m|a-1) 

Corresponding current will not depend on {Pb—Pa) since 
this part of the flux is removed from the LF. The current 
is therefore obtained by multiplying the LF by density 
of steps. Pa and Pb are relative probabilities so that 
Pa + Pb = 1- Accounting for this effect the expression 
for the current becomes. 



n\m\F{PB ~ Pa) 
2(l + |m|)(;c"i + |m|a-i) 
n\m\F 

~ {1 + \m\){lc' + \m\a-^) 
\m\ 

2{1 + \m\){l-^ + \m\a-^) 



In the limit of small m, the current reduces to. 



j(x) = h{PB - PA)\m\lj2 - hl-^dx{m^F) 



(3) 



(4) 



The second term in general will have the form V(Vft.)^. 
This term was derived using the Burton-Cabrera-Frank 
(BCF) theory and assuming that at small slopes the par- 
ticle density on the terraces depends on the even powers 
of local gradient (|J^. In the presnt case the difference 
in the velocities of two consecutive steps on a slope is 
related to the difference in the adatom densities on their 
terraces. However we do not explicitly assume the depen- 
dence of density on the even powers of slope. Further, in 
ref. Q it was conjectured that such a term can arise due 
to the differences in the velocities of the steps near the 
top and the bottom of a profile. 



In the hmit of large slope, \m\a 
reduces to, 



» L the current 



PA)aF fiFa^ 
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(5) 



The second term is proportional to -^dxTn, which was 
derived in the large slope limit in reference |^ . Thus the 
geometrical dependence of the symmetry breaking term 
in Eq.(3) exactly matches with previously derived two 
terms in the small and large slope limits. This shows 
that present method of deriving current, appropriately 
accounts for the physical processes contributing to it. 

We further argue that a curvature dependent cur- 
rent must be present in any adatom relaxation pro- 
cess that involves downward hops across the descending 
step edges. This argument is based on the observation 
that, in a 1+1 dimensional simulation, if adatoms are 
restricted completely to the in-plane hops (infinite SE 
barrier) then correlations do not grow beyond the dif- 
fusion length. On the other hand, when such hops are 
allowed, correlation length for stable growth and mound 
size for unstable growth increases in time 1 1 1 1 . Hence if 
the downward hops are allowed but the current is tilt- 
independent, then it may be expressed as a linear combi- 
nation ai'V^h + a2\'^h + ... including the nonlinear terms 
of the form V(V'^/i)^. We will retain only V'^h in the 
current corresponding to our minimal growth equation. 
Thus the form of the current corresponding to the mini- 
mal growth equation is 



j(^) = 



n\m\F{PB - Pa) 
2(l + |m|)(C' + |m|a-i) 
hlmlF 



dx 



2{1 + \m\){lc^ + \m\a-^) 
ItoI 



{l + \m\){lc' + \m\a-^) 



(6) 



The first term has been studied widely as the stable 
growth mode ||l^,|l^ and as the unstable growth mode 
Jl^ . In order to find the effect of asymmetry term on 
the growth exclusively, we have performed simulations 
of a 1-1-1 dimensional solid-on-solid model with no diffu- 
sion bias. A fourth ordered equation was earlier proposed 
by Villain |^ for similar situation. Under this condition 
Pa = Pb and first term vanishes in the Eq(3). The resul- 
tant growth equation in the moving frame with growth 
front is of the form, 



dfh = vV^h 



to 



{l + \m\){{le' + \m\a-^) 



(7) 



where, i^a is the appropriate constant for the asymme- 
try term and rj is white noise with < r]{x' ,t')r]{x,t) >= 
D5{x' — x)S{t' — t). In the limit of small slopes, renor- 
malization group (RG) analysis shows that the rough- 
ness exponent a = 1 and the roughness evolves with 
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the exponent (3 = 1/3 Q. For large slopes, asymmet- 
ric term, on power counting leads to z — 3a, while noise 
term gives z — 1 — 2a. If i/q is assumed to be invarient 
under renormalization, then again a = 1 and /3 = 1/3 
as in the case of small slope. Corresponding model is as 
follows. Atoms are rained on a 1-d substrate of length 
L randomly with constant flux. On deposition a given 
adatom is allowed to hop n times, as in a random walk. 
If the hopping adatom encounters another adatom be- 
fore n hops are exausted, the adatom stays there perma- 
nently. If n hops are exausted without any encounter, 
it stays permanently at the last position occupied after 
n hops. Thus for hopping, all sites with zero neighbor 
are equivalent, ensuring Pa = Pb- We have measured 
width W2, height- height correlations G'(r,t) and skew- 
ness a for this model where, W2 = Tlii^i — h) t^^ 
andG(r,i) = ^ X;r'('»(i"+r', ^))^- The skewness 

cr = w^/wy^ , where W3 = ^ Y^ii^i -hf H]. 

Presence of V^h term is verified from the flatness of 
the saturated width for small L. We have chosen, n=25 
giving Zc ~ 5. The saturated width is flat almost up 
to 5/c, showing that V^h dominates at small lengths 
Fig. ||(a) shows the morphology of the interface after 
80000ML are grown. As predicted by the Eq.(7), the 
asymmetry is evident in the flgure with a = — 0.47±0.05. 
. Fig. ^b) shows plot of W2 in time. We obtain initially 
/3 around 0.3 that attains a steady value of 0.355 ± 0.015. 
Initial small value can be related to the small slope re- 
gion that predicts a value of 1/3 (compare data in the 
region from 10 ML to 200 ML in the flgure with the line 
having slope of 2/3). Correspondingly h-h correlations 
lead to the roughness exponent a that increases from 0.5 
to 0.65 ± 0.01. These results indicate that most of the 
morphological features of the growth with diffusion with- 
out detachment are captured by the current expression in 
Eq. (6) . It also shows that diffusion of the adatoms rough- 
ens the growing surface. Diffusion bias causes additional 
effects in terms of stability or instability of growth. In 
particular, if the bias is varied from extreme -ve SE bar- 
rier to extreme -l-ve SE barrier, a stable— > unstable tran- 
sition is observed. In this transition however h — > —h 
symmetry is broken asymptotically. Note that for -ve 
SE barrier V2S/'^h term dominates with -|-ve value of 1^2 
[ p2t , so that asymptotically, asymmetric term becomes 
irrelevant rendering tr = 0. At exactly zero SE barrier, 
finite value of cr(« —0.5) is obtained, a can be regarded 
as the symmetry parameter, that changes abruptly at 
the transition point. Thus the growth transition is like 
2nd order phase transition. Details of such a transition in 
computer simulations of suitable models will be discussed 
elsewhere. 

Extension of this equation to 2-1-1- dimensions is possi- 
ble by similar kinetic considerations. For isotropic diffu- 
sion, same form as Eq.(6) is obtained with n = Vh/\Vh\ 
and replacing length-derivative product in asymmetric 



term by 



{l^^+\\7h\)\Vh\ ' 



For small slopes one obtains 



V(V/i)^ term which seems to describe many experi- 
mentally observed growth roughness measurements from 
vapor [ |l6[ . 

As mentioned earlier, the present method for obtaining 
current from kinetic considerations appropriately brings 
out the geometrical dependence in growth equation. We 
have applied this method to one of the stochastic growth 
models proposed to capture the essential feature of low 
temperature molecular beam epitaxy (MBE) the DT 
model 10|. Based on noise reduction technique, the 
simulations of this model |l7|| confirm that, 1) exponent 
(3 — 3/8, 2) the morphology is asymmetric with a « —0.5 
and 3) the current is tilt independent |§,|l3. The ob- 
served (3 value of 3/8 is constant almost over 8 orders of 
magnitude [T^ ]. The relaxation rules for adatom in this 
model allows it to hop only when it is deposited at site 
A or B (see Fig. |l|). Also only downward hop is allowed, 
if deposited at A and towards the step, if deposited at B. 
If it has choice of sites A, A or B,B or A,B etc. on two 
neighboring sites, then it will hop randomly to the left or 
right. Applying the considerations for obtaining current 
for this situation shows that the LF approaching sites A 
or B in the present set of rules is ,f , -1 , since the 
flux is affected only at large enough slope when site A 
and B overlap due to the short terrace width. Further, 
the LF is affected by relative motion of steps only when 
sites A and B differ by a lattice constant. The expression 
for the current is then. 



Jdt{x) = 



d'^m fiFa? 



m 



dx 



(l + |m|)2 



(8) 



From the rules it is clear that P4 = Pb- The growth 
equation corresponding to this current in the moving 
frame will be. 



l + \m\ (1 -f |m| 



+ 7? 



(9) 



where, Ua accounts for various constants in the corre- 
sponding expression for current. The power counting in 
this equation leads to , 2; = 4 from the first term, and 
z — 1 + 2a for large slopes corresponding to the second 
term which is expected to be operative mainly over large 
inclinations. The relation obtained from the second term 
is exactly same as the one obtainable from the noise term 
77. For z = 4 all the terms are marginal. This implies that 
z = 4 and /3 = 3/8. The second term is symmetry break- 
ing. Hence, above equation accounts for all the observed 
facts mentioned above in the simulation of DT model. In 
2+1-dimensions, with above rules for adatom relaxation, 
the local density of sites A and B need not be equal since, 
fluctuations in step edges render configurations that show 
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bias for sites A or B. As a result, slope dependent cur- 
rent will dominate the growth [T7| . Further implications 
of this equation are presently under study. 

In conclusion, we have proposed a simple method for 
obtaining current in a solid-on-solid growth in 1-1-1- di- 
mension. The resultant equation shows that presence of 
diffusion alone is resposible for roughning of a singular 
surface. It induces an asymmetric term in the contin- 
uum equation. As a result a stable to unstable transition 
will be associated with symmetry breaking. We have also 
demonstrated the application of this method of obtain- 
ing current to a stochastic model referred as DT model. 
The corresponding continuum equation is shown to pre- 
dict observed exponents and the nature of the interface 
profile correctly. 
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FIG. 1. A typical step structure formed during growth 
along positive slope, v and v' are velocities of the steps . 
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FIG. 2. (a) Morphology of the surface after 80000 number 
of layers. (b)Plot of width as a function of time. Straight 
lines with slope 3/4 and 2/3 are drawn for reference. 
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